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ABSTRACT 


A  model  for  the  interface  of  two  media  with  different  background 
densities  ,  separated  by  a  charged  bilayer, is  solved  exactly 
in  a  two  dimensional,  one  component  plasma  at  reduced  temperature  2. 

when  the  two  media  are  in  direct  contact  (no  gap) ,  then  we  can  think 
of  it  as  a  model  of  a  classical  metal  and  an  electrolytic  solution 
or  just  two  metals.  The  contact  ootential,  or  potential  of  zero  charge 
apoears  to  be  a  very  simple  function  of  the  ratio  of  the  background  densities 
We  also  find  that  the  potential  does  not  depend  on  the  surface  charge, a 
fact  that  we  exolain  using  a  corollary  of  the  perfect  screening  theorem. 

A  second  case  of  interest  is  tha  case  of  the  two  media  separated  by 
a  gap  that  in  our  case  could  be  charged  :this  is  a  model  of  a  two  dimensional 
membrane, or  the  oolarization  layer  of  the  inner  Helmholtz  region. 

Very  surprisingly, we  find  that  also  the  surface  dipole  is  completely 
screened  by  the  charges  surrounding  the  interface, so  that  the 
potential  drop  across  the  interface  only  depends  on  the  logarithm 
of  the  charge  densities,  as  found  bv  Ballone , Senator®  and  Tosi 
for  the  most  simple  case  of  a  discharged  contact. 


I -INTRODUCTION 


One  of  the  most  interesting  exact  results  in  the  theory  of  charged 
interfaces  in  recent  times  is  the  solution  of  the  one  component  plasma  (OCP) 
in  two  dimensions  by  Jancovici  (1,2).  Although  this  model  is  exactly 
solvable  only  at  reduced  temperature  2, which  in  a  three  dimensional  world 
corresponds  to  a  little  too  concentrated  ionic  solution  (for  biological 
systems), it  may  be  solved  for  a  rather  wide  variety  of  inhomogeneous  systems 
One  of  these  is  the  interface  between  two  media  of  different  background 
density, separated  by  a  charged  gap.  This  case  could  be  a  model  for 
the  classical  metal-metal  junction, the  semiconductor  junction 
the  metal  electrolyte  interface  among  others. 

In  Biology  there  is  also  a  system  which  corresponds  to  this  model: 

The  Donnan  equilibrium  is  established  between  two  media  containing 
different  concentra -ions  of  oroteins,  and  which  are  seoarated  by 
a  membrane  that  allows  passage  of  the  small  ions  only  (3) .A  system  of 
Darticular  interest  is  the  nerve  membrane (4 ,5) .The  mechanism  of  oroduction 
of  the  so  called  action  potential  is  strongly  related  to  the  charge 
distribution  across  the  inter  faces. Because  of  the  small  size  of  the  system 
and  the  low  conductivity, dir ect  electrochemical  measurements  are  difficult 
Only  recently  experiments  involving  interacting  double  layers  in  a  system 
similar  to  a  Donnan  system  have  been  reported  (o).iere  two  immiscible 
electrolytes  are  studied  electrochemical ly.dowever  these  experiments 
involve  rather  complex  organic  molecules,  and  it  is  very  nard 
to  construct  a  microscooical iy  correct  tneory  for  the  charme  anm 
and  ootential  profiles, and  tne  differential  caoacitances.  »e  should 
mention  in  this  context  the  recent  wcrx  of  Levine  and  Cuthwaite  (7) , 

Jicncils  anc  Pratt  (o)  and  Allastuey  and  Levesaue  (9). 


The  ouroose  of  this  work  is  to  oresent  a  simple  model  for  these 


systems  based  on  the  two  dimensional  OCP .For rester  and  3mith(JLG) 

have  solved  the  case  of  a  plasma  contained  between  two  fixed  walls, and  Smith 

solved  the  case  of  a  flat  ideal  electrode  with  image  forces  (11) 

Our  method  of  solution  is  based  on  the  method  of  Jancovici  (1,2) ,but 
follows  the  oultline  of  these  last  two  references. 

It  is  our  hope  that  this  exact  solution  will  serve  as  a  bencnmark 
to  evaluate  the  accuracy  of  different  aproxim3te  theories, and  to  check 
the  validity  of  exact  theorems  such  as  the  perfect  screening  sum  rules. 

In  section  2  we  give  a  description  of  the  model  and  a  brief  outline 
of  the  method  of  solution. The  reader  interested  in  the  technical  details 
of  the  solution  should  consult  references  (1,2,10,11).  In  section  3 
we  discuss  the  case  in  which  tne  width  of  the  gap  is  zero.  In 
section  4  we  present  calculations  for  the  case  in  which  the 
memorane,  or  gao,  is  of  finite  width.. Se  remark  that  this  case 
also  represents  the  inner  Helmholtz  layer  with  a  fixed  dipole. 


i 


1 


1 1 —■’IE T HOD  OF  SOLUTION 

The  system  consists  of  an  empty  strip  of  width  £,  that  separates 

the  two  plasmas  .  The  neutralizing  background  in  each  side  is  of  density 

-  eoc  /O'  and  —  eoe/fT  ,  where  e  is  the  elementary  charge,  and 
1  2 

the  borders  of  the  strip  have  charge  densities  -eo//1"T  and-ec'/TT 

1  2 

The  plasma  coupling  parameter  is 

where  p  =  l/kT  is  the  Boltzmann  thermal  factor, T  being  the  absolute 
temperature.  The  value  I"1  =2  may  be  special  in  the  sense  that 
the  pair  correlations  of  the  bulk  have  gaussian  rather  than  exponential 
screening .Otherw ise  this  parameter  is  not  unphysical , and  , in  fact 
it  corresponds  to  the  order  of  magnitude  of  a  1.1  electrolytic 
solution  at  room  temperature. 

Following  Jancovici  (1,2)  we  consider  initially  our  system  to  be 

confined  to  a  disk  of  radius  H.  In  this  disk  there  is  a  ring  of  inner 

radius  F  and  outer  radius  F  .  The  uniform  oackground  density  of 
1  2 

the  inner  region  is  -eaC/TT,and  that  of  the  outer  region  is-eCC/TT 

1  < 

Clearly 

a  +£  (2.1) 

2  1 

where  £  is  the  width  of  the  ring  which  eventually  will  become  our 
membrane. In  our  model  this  memebrane  will  be  just  a  mao.  The  mooile  ions 
crane  e  are  free  to  be  anvwhere. 


The  total  number  of  ions  .V  must  satisfy  the  electroneutrality  relation 


2  2  2 

N=2(tf  R  R  )+R  CC+(R  -R  )  oc  (2.2) 

1122  11  12  2 


The  hamiltonian  of  this  model  is 
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2  )  2  ")  2 

H»  (1/2)  e  {  -  lnlr  -r  (  +  [  2or"  R  ln(r  /R  )  d(r  -R  ) 

k>j  k  j  k=l  1  1  k  1  k  1 


2  2  2 

2tTP.  In  (r  /R  )  S(r  -R  )  +  oc  ( r  -R  )  d(R  -r  ) 
22  k  2  k  2  1  k  1  Ik 

2  2 

+  oCR  In  (r  /R  )  e(r  -R  ) 
11  k  1  k  1 

2  2  2  2 
+  oC(r  -R  )-acR  In  (r  /R  )  0(r  -R  )  )  }  +  3 

2  k  2  2  2  k  2  k2  H 


(2.3) 


where  r  is  the  oosition  of  ion  k,  of  charge  e,and  w(x)  is  the  Heaviside 
k 

function  of  x.B  is  a  background  term  that  is  irrelevant  to  our  present 

M 

calculation .The  potentials  have  been  shifted  to  insure  continuity 
across  the  olates. 

<ve  comoute  the  canonical  oartition  function 


Z  = ( 1/M ! 

:i 


Replacing  (2.3)  into  (2.4) 


(  2.  4) 


-e~;. 

(2TT) 


z  r  +1 
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Z  £. 

-QC  (r  -R  ) 
i  1 


2(n  —  n)  R  2n  +1 

9  (  2  g 

i-  [1/  R  ]  dr  r 

1  Jr 
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2  2 

2 ( n  -n)  2(n  -n  )  R  2n  +1  -of  (r  -R  ) 

g  2  g  f  2  2  2 

*  [ 1/R  ] [1/  R  ]  dr  r  e 

1  2  Jr 
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where 


(2.5) 


n  »n  -2<f  R  -  ct  R 
g  1111 


(2.6) 


2  2 

n  =  n  -2<T  R  -2cr  R  +(  oc  R  -  oC  R  ) 
2  11  22  22  11 


(2.7) 


and  p=(l/kT)  is  the  Boltzmann  factor. 

Introducing  the  incomplete  gamma  function 


i,b)=  J 

"'a 


b  a- 1  -  t 

dt  t  e 


(2.8) 


we  qet 


.>1  VI 

(  ft) 
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11  'll  1 
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g 


2 ( n  +1) 

g 


+  [1/  R  ]  [  1/2 (n  +1)  HR 

1  g  2 


-R 


2(n  -n)  2(n  -n  )  2 

9  2  g  2 

(1/R  ]  [  1/  R  ]  [  o  ( n  +1J0CR  )  -tf(n  +1;  oc 

1  2  2  2  2  2 


n  +1 

2  2 

expfccR  ]  (1/at  )  } 
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2  2 


The  single  particle  density  is 


N- 1  2  2  2  n 

(1/D  )  {  exp[-ac  (r  -R  )](r/R  )  9(  R  -r 

n  1111 

n*0 


2  n 


(r/R  )  [©(  r  -R  )  -©(  r  -R  )  )  + 

1  12 


2  n  2  n 

2  g  2  2 

+  (r/R  )  (R  /R  )  exD(-oC(r  -R  )  ]  [9-(  r  -R  )  ]  } 

2  2  1  2  2  2 

(2. 


with 


2  2n  2 

D  =exo  (OCR  )  (1/R  )  #(n  +l;acR  )  + 

nil  1 


2 

+  (v'  /  ( 2n  +2) 
i  q 


2  n  +2 
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( R  /R  ) 
2  1 
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2 

R  )]  I 
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2 ( n  -n  ) 

2  g  2 

(1/  R  )  [  K(n  +l;Srf) 

2  2  l 


-&{r\  +1;  R  #)} 
2  2  ■*- 


n  +1 
2  2 
exp[ac  R  ]  /at 
2  2  2 


(2.11) 


Similarly, the  pair  density  distribution  function 


is 


j  m  m  m  m 

p(r  ,rj=  p(r  )f(r  )-/_ _ ,  r^r^  \  1  (A 


12  '  1  2  m  ,  m 

1  2 


2  (  m  +m  ) 
2  12 
[1/R 
2  V  2}  1 


{  G  (  r  )  3  (  r  )  }  /  { D  D  } 
1  2  mm 

1  2 


(2.12) 


where  D  is  defined  in  (2.11)  and  3(r)  is  given  bv 


2  2 

:  ( r)  =  exD[-oc  (r  -R  )]  e(  R  -r  ) 
1  1  1 


)  {[9(r-R  )-e(r-R  ) ] 

1  1  2 


‘A  1 A 

q  *"*2  ^ 

*  (R  /;?  J  )  exo  [  -oc  ( r  -R  ))  -r) 

2  1  2  2  2  2 


In  the  limit 


R,R  fR  ->  oo 
1  2 

R  -R  (2.14) 

2  1 


and  using 


^(n  +1;N)  =  (V2TT  n/2)  exp(-r*n  In  n)  {1  +  ^[(N-nl/Jin  ] 
with 


(2.15) 


2 

(1/D  )  {exo(-2cf  x  -2xtfT  )  9(-x)  + 

t  1 


exc  (-2x<(2‘(  t+  la'll  [  0(x)-0(x-E)} 
1 


+  exo[-2ae  (x-  t)  -2  (x- C)  >/2(t+\j2cr  )-2$(2cf  +£t)  ]  9(x-t)J  } 


with 


D  *exp(  t  /oc  )  [  1  +<F(t  /Joe)) /{of 
t  1  11 


-  [exp{-2ff?(t+v20*  )}-l]/[FT(t+  {ler  )] 

1  1 


+  exp{-2f  i|7(  tW^cf  )  +[  t+>I7(o* +cr  )  ]  /oc  }  [  l-i{  t+  (7(cr  +  p  )/ Jor  } ]  /{ac 

1  122-*-  1222 

(2.17) 

and  for  the  pair  density 

2  2 

c(z  ,z  )=  p(x  )p(x  )  -exp[-2cc  x  -2cx  x  ) 

J  1  2  '  1  2  i  j 


3  roo  2 

4/ft  I  dt  exp[-  2t(x  +x  +  -iy)  ]  /D  I 
J -oo  12  t 

(2. Id) 


where  z=x  +  iy, a  result  due  to  Jancovici  (private 


communication, 12) . 


Ill  THE  3ACKGR0UND  JUMP  MODEL 


A  very  interesting  particular  case  of  (2.16)  arises  when  the 
width  of  the  interface  layer  £  is  zero.  In  that  case  the  charge 
density  is 

( _  rOO  2 

p(x)=  (2/^fT)  dt  (1/D  )  {exp(-2ccx  -2xt \[2  )  x<0 

J  J- 00  t  1 


o(x)  = 


(1/0  ) 
t 


2 

exp(-2ocx  -2x\f?(  t+\/7cr)  ] 
2 


x>0 

(3.1) 


D  =  exp(  t  /  cc  )  (  1  +  |  ( t//cc  )  ]  /^cc 


+exp((t+^cf)  /  cc.)  [  1  -  l’(t/N/oc  )  ]  //oc 

2  x  2  2 

(3.2) 

This  equation  for  the  case  <3=0  has  been  independently  obtained 

by  Jancovici  (  to  be  published  ,12). 

Figures  2  shows  the  mobile  charge  orofile  for  various  ratios  of 

the  backcrround  densities  cc  /cc  .The  value  15  is  scaled  to  reocesent 

2  1 

the  difference  in  density  between  a  1  molar  electrolytic  solution 
and  a  simple  metal:  If  we  take  this  number  to  be  64,  then  to  go 
from  three  dimensions  to  two  dimension,  we  simply  take  the  (2/3)  oowe 


which  is  16. Clearly,  most  metals  are  far  from  being  classical 
one  component  plasmas  (Drude  Theory) , although  semiclassical  treatments 
are  widely  used  to  explain  transport  properties  (14) , which, as  we 
know,  are  strongly  dependent  on  cor  relations. Unfortunately  the  quantum 
mechanics  of  fully  correlated  systems  is  very  difficult, more  so  in 
the  vicinity  of  a  surface. For  this  reason  we  believe  that  an  exactly 
solvable  model  of  a  fully  classical  metal  interface  is  useful  because 
it  illustrates  the  subtleties  of  the  behaviour  of  correlations  in 
charged  systems. Figure  3  shows  the  total  charge  distribution  near 
the  interface  as  a  function  of  the  background  ratio. Although  it  is  clear 
that  the  penetration  depth  on  the  metal  side  is  smaller  for  the 
'metal'  side  (  side  2), the  charge  densisty  is  also  much  larger 
Figure  4  and  5  show  the  rather  strong  dependence  of  the  charge 
profile  on  the  surface  density  6" .  In  Figures  6  and  7  we  have 
compared  the  charge  profiles  for  different  O',  but  keeping  o' constant. 
Again,  we  see  the  rather  dramatic  dependence  of  the  charge  profile  on  the 
surface  charge.  Figure  ti  shows  the  dependence  of  the  contact  potential 
on  the  surface  charge  O', compared  to  the  ideal  electrode  case. 

To  investigate  further  this  point  consider  Poisson's 
eauation  (15) . 


where  oc(x)=oc  for  x<u  and  ac(x)=ccfor  x>0. 
1  2 

Sy  elementary  integration 
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( x ) I x I dx 
'  -oo 


(X  -x)  [p(x  )-oc(x  )  J  -TT 

1  J  1  w 


(3.4) 


Using  (3.1),  we  get  after  some  straightforward  calculation 

J(oo)  -(j>(-oo)  =.(1/4)  ln(QC/0C)  (3.5) 

This  extremely  simple  formula  for  the  contact  potential  of  two 
metals  (or,  cor  respond igly,  to  the  potential  of  zero  charge  for  the 
case  of  the  metal  electrolyte  interface) ,  is  not  really  unexpected 
when  the  charge  o=0:Actually  this  result  was  obtained  on  the  basis  of  a 

2 

f  =F  /kTN=—  (1/4)  ( e  /kT)  In  cc  +F  (T)  (3.6) 

i  i  i  o 

where  i  is  the  index  for  the  side  of  the  interface  (i=l,2). 

The  reversible  work  to  transfer  one  mobile  charge  from  side  to  the  other 

is  clearly 

2 

Jf=(-l/4)  (e  /kT)  In (ct  /oc  )  (3.7) 

2  1 

But  the  reversible  work  is  only  electrostatic, so  that 
2  2 

Af=(e  /kT)£lC  *  (-1/4)  (e  /kT)  ln(oc/oe)  (3. a) 

J  2  1 

where  we  must  remember  that  because  of  our  special  choice  of  units 
the  electrostatic  potential  is  measured  in  units  of  2kT^ 

This  result  is  contained  in  the  work  of  Ballons,  Senators  and  Tosi 
(la)  who  showed  for  the  neutral  interface  that  the  potential  drop 
across  the  interface  depends  only  on  the  difference  of  the  chemical 


Dctential  in  the  bulk  ohases. 


The  surprise  is  that  the  surface  charge  has  no  effect  whatsoever 
cn  the  potential  drop  across  the  surface. Clearly  this  must  be 
due  to  the  fact  that  because  of  the  perfect  screening  theorems  (16,19,20) 
the  charge  distribution  has  no  multipole  moment. Indeed  the  theorems 
have  been  proven  for  finite  size  disks  (  for  R  finite  ),and  in  that 
case  <x>=0, which  implies  no  potential  drop  due  to  charge  &  .we 
conjecture  that  this  relation  remains  valid  in  the  limit  or  R->oo  . 

Just  as  an  interesting  aside,  if  we  scale  the  relation  (3. b) 
to  three  dimensions, then,  since 

(2/3) 

f>  -n  (3.  9) 

A(0*  (-1/6)  ln(  n  /n  )  +  F  (T )  (3.10) 

J  2  1 

where  n  ,n  are  the  bulk  electron  densities  of  the  metals, and  F(T)  is 
1  2 


is  a  function  of  the  temperature  alone. 


IV  THE  GAP  MODEL 


The  equation  (2.16)  offers  a  rich  variety  of  possibilities  that 

we  shall  not  explore  at  this  time. We  have  plotted  the  charge  density 

densities  for  two  different  cases  in  figures  9 

and  (10)  (total  charge  density) .More  complicated  profiles  are  obtained 

when  the  interface  is  charged  and/or  polarized.  However,  the  more 

interesting  quantity  is  the  electrostatic  potential , which  is  obtained 

by  integrating  Poisson's  equation  (3.3).  The  results  for  the  charged  gap 

where  cr  =o*  =o'  and  the  background  <z  =4  are  shown  in  Figures  11  and  12 
12  2 

again  here  we  observe  that  the  effect  of  the  charges  is  screened  out 

completely,  and  the  potential  drop  is  given  by  (3.5).  We  finally 

computed  the  case  in  which  the  interface  had  a  permanent  dipole 

by  setting  <r'=-o'=2.  Figure  13  shows  the  two  oonosinq  orientations 
1  2 

of  the  dipole,  just  to  show  that  they  are  also  completely  screened  out 

Figure  14  is  the  same  calculation  but  with  qc  *4. Again  here,  the 

2 

potential  droo  is  given  by  (3.  5).  This  is  a  consecuence  of  the 
perfect  screening  theorems. 
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FIGURE  CAPTIONS 

Figure  1  System  geometry. 

Figure  2  Mobile  charge  density  profile  fk)(x)  Curve  1,  oc  =1;2  ,  oc -2;  3  ,  ot 

3  2  2  2 

4,OC  =  d;5,OC  =  lo 
2  2 

Figure  3  Total  charge  profile  Tf  p(x  )-odx)  ,same  numeration  of  curves  as 

Figure  2. No  gaD  in  this  case*f=Q 

/ 

Figure  4  same  as  Figure  3  but  with  d* 1.  . 

Figure  5  same  as  Figure  4  but  with  d*-l. 

Figure  6  Variation  of  total  density  orofile  with  charge  for  cc  =  2:Curves 

2 

1,2, 3,4, 5  correspond  to  surface  charge  densities  0^-4  ,-2  ,  J  ,2 , 4 

Figure  7  same  as  Figure  o  but  with  or-=16 

2 

figure  a  Contact  density  as  a  function  of  the  surface  charge 

Curve  2,ac  =4;  3  ,cc  =  lo 
2  2 

Figure  9  Total  charge  density  profile  for  a  gap  model  with  £=1 

Here  noth  elates  are  equally  charged  Curve  1 ,d  =2;Curve  2  no  charge 

i  2 

Curve  3,  d*-2. 

Figure  ID  3ame  as  figure  9,  but  with  O'  =-6  =  +  2  , oc  =4 . Case  with  dioole 

12  2 

and  background  gap. 

Figure  11  Potential  profiles  for  the  same  cases  as  in  Figure  9. 

Fiaure  12  Same  3S  figure  11,  but  with  a  =4. 

2 

Fiiure  13  Interface  with  dioole, no  background  gap. Potential  for 

same  case  as  Figure  1J  with  cc=l. 

2 


Figure  14  Same  as  Figure  13, but  with  background  gap  cc =4 .Cor resDonds 

2 


to  the  charge  density  of  Figure  10. 
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